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Abstract 

In this work we study the spreading effects of wave packets on the two-flavor neutrino oscillation probabilities in vacuum 

for describing neutrinos by Majorana fields. To do it, the flavor neutrinos are considered as superpositions of neutrino mass 

r^ ■ eigenstates which are described by localized wave packets. The spreading effects are studied in two cases: (i) Considering only 

^vj ' the space spreading; (ii) considering simultaneously the time and space spreading. For the space spreading case, we obtain a 

, functional form of the oscillation probabilities which is consistent with that reported in the literature for Dirac neutrinos. For 

f) I the time-space spreading case, we obtain a new expression for the oscillation probabilities that includes some small changes in 

^^ ' the oscillation and coherence lengths. 
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I. INTRODUCTION 



Knowing the nature of neutrino fields is an open problem in particle physics [lj-|3| . This problem might be solved, 
experimentally, establishing if neutrinos are: (i) Majorana fermions; (ii) Dirac fermions. In the first case, as a 
consequence, neutrino and anti-neutrino are the same particle being described by two-component spinorial fields 
called Majorana fields [i|-y]. In the second case, neutrinos and anti- neutrinos are described by different four-component 
spinorial fields called Dirac fields [iHS]- However, independently of its nature, the theoretical description of neutrino 
oscillations should lead to the same results for the cases in which neutrinos are described by Majorana or Dirac 
fields. As an example of this fact, for the case in which neutrinos were described by Majorana fields, we obtained the 
standard plane wave expressions for the two-fiavor neutrino oscillation probabilities [J] , which have been also obtained 
describing neutrinos by Dirac fields. 

The standard plane wave expressions for neutrino oscillation probabilities have been obtained in the context of 
several quantum mechanics treatments (for instance, see (5|-[l0j). On the other hand, in the context of different 
quantum field theory treatments, neutrino oscillations in vacuum have been extensively studied describing neutrinos 
by Dirac fields |lll - l23l |. For Dirac neutrinos, the space spreading effects of the mass eigenstates wave packets on 
the neutrino oscillation probabilities have been widely investigated |ll| - |2]| . Additionally, for Dirac neutrinos, some 
aspects of the time-space spreading effects of the mass eigenstates wave packets were studied in [IJ, [T^l ■ Until now, 
however, the spreading effects of the mass eigenstates of the wave packets on the neutrino oscillation probabilities 
have not been studied for neutrinos described by Majorana fields. 

The main goal of this work is to study the spreading effects of wave packets on two-flavor neutrino oscillations 
in vacuum for considering that neutrinos arc described by Majorana fields. To do it, we calculate the two-fiavor 
neutrino oscillation probabilities in the context of a wave packet extension of the quantum field theory treatment 
that we previously developed for the case in which the Majorana neutrino mass eigenstate were described by plane 
waves [^. In the present treatment, the fiavor neutrinos are considered as superpositions of neutrino mass eigenstates 
described by localized wave packets. Thus, the spreading effects of the mass eigenstate wave packets on the two-flavor 
neutrino oscillation probabilities are studied in two cases: (i) Considering only the space spreading; (ii) considering 
simultaneously the time and space spreading. For the first case, in the relativistic limit, we find that the Majorana 
neutrino oscillation probabilities have the same functional form as those reported in the literature for Dirac neutrinos 
[lll42l| . On the other hand, for the second case, in the relativistic limit, and working in the regime of oscillations with 
transversal spreading, wc obtain a new expression for the two-flavor neutrino oscillation probabilities that includes 
some small changes in the oscillation and coherence lengths. Specifically, we observe for this second case that the 
oscillation and coherence lengths are a little smaller than the ones obtained in the space spreading case. 

The content of this work has been organized as follows: In section two, we present for neutrinos described by 
Majorana fields how is possible to calculate the two-fiavor neutrino oscillation probabilities in the context of a quantum 
field theory treatment, for which the mass eigenstate are described by plane waves; in section three, we extend the 
plane wave treatment presented in the previous section, to the case in which the Majorana neutrino mass eigenstates 
are described by localized wave packets; in section four, we study the space spreading effects of the mass eigenstate 
wave packets on the two-flavor Majorana neutrino oscillation probabilities; in section five, we study the time-space 
spreading effects of the wave packets on the Majorana neutrino oscillations; finally in section six we present some 
conclusions. 



II. NEUTRINO OSCILLATIONS USING PLANE WAVES 

The standard plane wave expressions for the oscillation probabilities between two flavor left-handed neutrinos were 
obtained in the context of a treatment developed in the canonical formalism of Quantum Field Theory for the case 
in which neutrinos were described by Majorana fields |J|. In this treatment, the flavor neutrinos were considered as 
superpositions of mass eigenstates with specific momenta and the mass eigenstates were described by plane waves. 
For the case of the relativistic limit (L ~ T) and after including a normalization constant, the standard plane wave 
expressions for the neutrino oscillation probabilities [J] are written as 



PJ'^(L) = l-sin2 [26*12] sin^ 
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where Amfj — "^2 ~ '^ii -^ is the energy of the neutrino, L is the distance between the neutrino source and the 
detector and sin^[20i2] is given by 



sin^[20i2] = 



4Ai 



(3) 



(1 + Ai)2' 
with 012 representing the mixing angle between the two mass eigenstates in the vacuum. In the last expression Al is 



where R^ is defined by means of 
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and the parameters ttj^^ , nii, ^ and tti^^ ^^ are related with the masses mi and m2 of the neutrino fields vi and V2 
(with definite masses) through the following relations 
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being ul the Major ana complex phase [J]. 

The expressions ([T]) and ([2]) are obtained starting from the oscillation probabilities defined by 
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in such a way that in the space-time production point (xq) the initial left-handed neutrino flavor state (electron 
neutrino flavor) i^f"'(a;o)) is deflned by the following superposition of the mass eigenstates |:^f"'(xo)) and \i^P'^{xo)) 
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where a sum over helicities is taken in the superposition. The mass eigenstates |i^f"'(a;o)) and \i'P^{xo)) involved in 
(1101) are obtained using plane waves from the vacuum state |0) as 



Wr{xo))=Ae^P^^»aiip,,h)\0) 



(11) 



where A is a normalization constant, aj is the creation operator of a Majorana neutrino of deflned mass, xq is the 
space-time point where this Majorana neutrino is created, Pa = {Ea,Pa) is the four-momentum of the mass eigenstates 
and a = 1,2. We have assumed that each mass eigenstate involved in ([TU|) has associate a specific four- momentum. 

In the oscillation probabilities ^ and ^ , the fiavor neutrino field operator i)^ is defined as a superposition of field 
operators of neutrinos with defined mass i>a by means of the expression 



i>a{x) = ^ UL^^i>a{x), 



(12) 



where a = cl, fJ-L and Ul is an unitarian rotation matrix given by 
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The field operators of Majorana neutrinos with defined mass Va involved in ([T^ are defined as [j| 
d^p 

h=±l 
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where E'^ = \p\ -\~ m'^ is the energy of the neutrino field with defined mass and x^{p) a-re Majorana spinors with 
helicity eigenvalues ±1. 



III. NEUTRINO OSCILLATIONS USING WAVE PACKETS 



In this section, we will extend the plane wave treatment for Majorana neutrino oscillations that we have have 
presented briefly in section two (details of the treatment can be seen in Q), for the case in which mass eigenstates are 
described by localized wave packets. In this treatment, we do not focus on the study of the details of the interaction 
processes in which neutrino is produced and detected. Here, on the other hand, it is assumed that wave packets 
describing mass eigenstates are localized and the coefficients of their superpositions are given by the elements of the 
unitarian rotation matrix Ul given by (|13p . The matrix Ul establishes a relationship between the flavor and mass 
eigenstates bases of Majorana neutrino fields in vacuum. 

There are different reasons for understand why the description of mass eigenstates using wave packets is most 
appropriate to study the neutrino oscillations with respect to the description from plane waves [l|-|3| . Some of these 
reasons are that the neutrino source and the detector are localized and there exists a spread for the neutrino momentum 
[l| . Additionally, we have to keep in mind that plane waves localized in some point Xq are in an obvious contradiction 
with the Heisenberg uncertainty principle. Here, we take into account these reasons when we describe the Majorana 
neutrino mass eigenstates in terms of superpositions of localized wave packets. To do the last, we first consider that in 



a point Xq 



(ig, rg) is created a left-handed electron neutrino which is described by the following superposition 



of Majorana neutrino mass eigenstates \v^''{xo)) and \u^^{xo)) 
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In contrast to the expression ([TT|) , now the mass eigenstates involved in P^ are written in term of localized wave 
packets in the form 



ii.r"(^o)> = A 
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where A is a normalization constant, ipa{p, (Pa)) is a probability density function which depends on the momentum 
Pa and the average momentum (pa), and a = 1,2. In general, ipaip,{Pa)) rnay take any form, but it is usually 
approximated by a Gaussian distribution assuming that is peaked around the average momentum [7|, Illl423l | . In this 
work, we define the probability density function as [J] 



Mp, (Pa)) - (27r)-3/4[Detr]i/4exp 
in such a way that ipa = i^a{p, (pa)) satisfies the conditions [3| 
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In the probability density function given by (J17p . we have taken the convention of summation over the Latin repeated 
index k and j. Additionally, we have assumed that the spreading over the mass eigenstates \h'^'^{xo)) and |i^"'^(a;g)) 
is the same, because these eigenstates are created simultaneously by the same weak production process. This fact is 
the reason that justifies why the matrix F^j is identic for both mass eigenstates. Moreover, it is important to note that 
this matrix is symmetric and the eigenvalues of its inverse FT. are the squares of the widths in the momentum space [3| . 

The oscillation probabilities between two flavor neutrinos using wave packets are 
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where i^Y^^ixo)) given by (fT5|) is a superposition of mass eigenstate wave packets \v^^^{xo)) and |j/f'^(a;g))) in the 
creation point. The expressions ([^0]) and (PT|) describe, respectively, the probabilities of finding an electron neutrino 



(i/e) and a muon neutrino [v^j) at a distance i in a time T of the creation point Xq. In the calculation of the transition 
probabilities (|20p and pip , the relativistic dispersion relation is approximated by means of an expansion around the 
average momentum of the wave packets {pa) [13 

Ea{p) ^Ea+Va- {p ~ (Pa)) + ^(p - (pa»'f^fc,(p - (pa))' + ' ' ' , (22) 
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Additionally, it is possible to write that 
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The highest power of {p— (pa)) in the expansion of the energy given by (|22l) determines the two different cases which 
we will studied below: (i) If the highest power is taken until second order, then the space spreading of the mass 
eigenstates wave packets will be considered in the propagation of the wave packets; (ii) If the highest power is taken 
until third order, then the time and space spreading of the wave packets will be simultaneously considered in their 
propagation. The time spreading effects on the neutrino oscillation probabilities are only observed in the second case, 
while the space spreading effects on the neutrino oscillations are observed in both cases. 

IV. SPACE SPREADING EFFECTS ON NEUTRINO OSCILLATIONS 

In this section, we will study the space spreading effects on the neutrino oscillations. To do it, we expand the energy 
given by p2p up to second order in the power series of (p— (pa)) [15|. This fact is justified in the fact that the width 
of the wave packets is very narrow around the average momentum. For this case, the matrix Tkj can be diagonalized 
by means of an orthogonal transformation, i. e. a rotation |3|, ll5| . Given that the expansion of the energy does not 
change this rotation, without lost of generality we can take a reference frame where the matrix is diagonal 

Tfe, = ^Sk,, (28) 

'^p 

with (Tp representing the width of the wave packets in the momentum space. We assume that the width has the same 
value for each of the dimensions of the momentum space, due to the wave packets are taken as isotropic. Additionally, 
we define the width of the wave packets in the coordinate space ar through the uncertainty relation 

(^rO-p = -. (29) 

If the energy given by (22) is substituted in (PH)) and (PT|) . keeping up to the second order in the power series of 
{p— (pa)), we obtain that the neutrino oscillation probabilities are written as 

P7-(r,L) = ^^J,y/, (^ ^\,^2 {A' cxphAi0f(r)] +exphA2</^f(T) + A^H exphA3C(r)]} , (30) 

P^^'-iT^L) = ^^J,y^, (^ ^^2)2 {^' exp[-Ai0f (T)] + A^ exp[-\,cl^UT) A^H exp[-A30f (T)]} , (31) 



where Ai = A2 = l/2cr^, A3 = l/4cr^, T ^ t — to, L ^ |r — ro| and the functions in the arguments of the exponentials 
are given by 

rAT)HL-viTr, (32) 

cj)t{T)^{L~V2Tf, (33) 

0f (T) = (i - vif + {L- V2Tf - lAa^Ei ~ E2)T + lAaHpi ~P2)L, (34) 

with Va = \va\- The quantity H that appears in the oseillation probabihties pop and pip is written as 

H ^ (^^ j ^l/2 E \/(-^i - M(Pi)l)(-^2 - fe |(P2)I). (35) 

Now, we take into account the fact that in the atmospheric and reactor neutrino oscillation experiments it is only 
possible to measure the distance between the neutrino source and the detector L, while the neutrino propagation time 
T is unknown [3, 113: [2l|- However, in the case of the accelerator neutrino experiments (for instance K2K, MINOS, 
OPERA) it is possible to measure the neutrino propagation time T [2l|l ■ By this reason, if we focus only on the case of 
atmospheric and reactor neutrino oscillation experiments, then it is necessary the elimination of the time dependence 
presents in ((30|) and ((3T|) . This last can be performed, if we take the average on the time of the expressions pop and 
([57)) in the following form 

P^r" (L) = / P.i'^" (T, L)dT, (36) 

P!:^'" {L) = j Pj-J^ (T, L)dT, (37) 

The time integrations can be performed using both Gaussian integration and the Laplace approximation method. 
After the time integrations are performed, we obtain from p6l) and p7p the following oscillation probabilities 



A4 



1/2 



^r^W=(rTA^ - + ^ + ^-1:^^) -M^f!-m. (38) 



1 Ta^ A2 / 2 \^^^ 1 

'(^) = 7T-TI^ - + - - 2A^S ^—3 exp [z/f - f^] ) , (39) 



where 



/f = {E, - E^f-^^L - (pi - P2)L, (40) 
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We have explicitly proved that if the average on the time in the expressions p6p and p7p is performed using Gaussian 
integration, the results are the same as the those obtained using the Laplace approximation method. In both cases, we 
have obtained the oscillation probabilities given by (|55)) and ([M]) . The functional form of the oscillation probabilities 
psp and p9p is in agreement with the one obtained in [7[, where it has been used a quantum mechanics treatment 
and the flavor neutrinos have been described by a superposition of mass eigenstates wave packets. 

In order to obtain from (|38l) and (p9| expressions for the oscillation probabilities in the relativistic limit, the 
following relativistic approximations are used [l^, ^M, 13 

E,c,E + ^^, (42) 

p^c^E+il-0^, (43) 

Va^l ^, 44) 

2E' ^ ^ 



where ^ is a dimcnsionless cocfBcient, typically of order unity, that depends of the neutrino production process and 
E is the neutrino energy determined by the kinematics of the production process for a massless neutrino. After the 
relativistic limit is taken, we obtain from (|38|) and p9|) the following neutrino oscillation probabilities 
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P7-(L) = isin2[2M|l-exp 
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where Lose is the oscillation length and Lcoh is the coherence length given by 
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in agreement with the corresponding lengths very well known in the literature [Tl, Illl42l| . To write the expressions 
(j45|) and ([46|) . we have used the definition of sin^[20i2] in terms of the parameter A given by ([3]), where ^12 is the 
mixing angle in the vacuum between the two mass eigenstates. Specifically, the probability (|45|) represents the survival 
probability that an electron neutrino (vg) be detected at a distance L in a time T of the creation point xq = (0,ro), 
where by simplicity to = 0. On the other hand, the probability ()46|) represents the probability of oscillation from 
an electron neutrino (i^e) created by the source at point xq to a muon neutrino (Vf^) measured by the detector at a 
distance L in time T. The dependence of the oscillation probability (|46p respect to the mixing angle 9l is in agreement 
with the reported in [23, [2^ , where the effects of positive and negative energy components of mass eigenstate wave 
packets on the two-flavor neutrino oscillation probabilities have been studied. 

The dependence of the oscillation probabilities P5)) and (P5|) respect to Lose and Lcoh is in agreement with the 
reported in the literature [3, [ll| - |2l| . The first term in the argument of the exponentials in (P5|) and pS)) is the 
standard oscillation phase proportional to the propagation distance L. The second term in the argument of the 
exponentials is called the damping factor [3] . This term implies a quadratical decrease of the oscillation probabilities 
with the distance L and determines how far the oscillations take place [l2l- For L ^ Lcoh, the interference of the 
neutrino mass eigenstates is suppressed and the oscillations due to Lose disappear J12l | . This behavior can be partially 
originated in the progressive separation of mass eigenstates wave packets propagating in space [l7| . The third term 
in the argument of the exponential is called the localization factor [l^l and this term implies that 
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meaning that the neutrino production process, which is characterized by the width of the wave packets tTr, is localized 
in a region much smaller that the oscillation length Lose [l2|. Thus, this factor does not depends on the distance L. 
As can be observed from (|T7)) and (|^, the oscillation length and the coherence length are related by 
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showing that the coherence length is much larger than the oscillation length |2l| . The maximum number of oscillations 
can be obtained from these lengths as 
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Because in the neutrino oscillation experiments one has L ~ Lose, then the term exp 
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is nearly equal 
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to one [21|. Additionally, it is easy to show that | wi — W2 | Lcoh 

is also nearly equal to one l2l|. In this form, the oscillation probabilities ([45]) and (j46|) can 
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that reduce to the standard plane wave expressions for the neutrino oscillation probabilities (H]) and ©. 

V. TIME AND SPACE SPREADING EFFECTS ON NEUTRINO OSCILLATIONS 

In this section, we will study the time-space spreading effects on the neutrino oscillations. This is possible if we 
expand the energy given by p2|) up to third order in the power series of {p— {pa})- For this case, we take a reference 
frame where the matrix T^j is diagonal and identical to (P5|) . Substituting the energy given by ^^ in ([20]) and (|2ip . 
keeping up to the third order in the power series of {p — (Pa)), wc obtain that the neutrino oscillation probabilities 
are written as 

K-iT,L) = j^j^ {^exp[-A.0r(T)] + ^exp[-A.C(T) + |^ exp[-A3C(T)] 

(54) 

K-iT,L) ^ (^ (TtW {^^^P[-^^^^(^)] + ^-phA.C(T) - |^exp[-A30nT)]} , 

(55) 
with H given by ([35]) . Ai = A2 = l/2cr^, A3 = 1/4ct^, and the functions in the arguments of the exponentials are 
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while the functions gf (T), g^iT) and gl{T) are 
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where the spreading time T^^'"' has been defined as T'^^p'' = -j-;^, with Sa = yf- The spreading time in neutrino 

oscillations was initially defined in the context of a quantum mechanics treatment in [7[. This time was posteriorly 

considered in the context of a quantum field theory treatment of neutrino oscillations in [ij, [l3| • 

We observe in ([M)) and (|55p that there are two different spreading times T^'"' and T^'"'. For simplicity, we will 
work in the limit in which the masses are nearly degenerate. In this case, it is possible to consider T/'"' = T^^"' and 
to work with only one spreading time defined by T^^^ ~ E^ /2a'^fh, with to the mass in the degenerate limit [17[. 
Additionally, we will work also in the regime of oscillations with transversal spreading defined for times T that satisfy 
Apa^/v < T < T,j,r [13]. If the transversal spreading is neglected (T^pr- ~ 00), this means that the time spreading 
effects are neglected in the oscillation probabilities, then the expressions ((54)) and (|55p lead coherently to (pO)) and 
(|5T|). respectively. 

For the limit of nearly degenerate masses and considering the regime of oscillations with transversal spreading, the 
oscillation probabilities ([M|) and ([55]) can be written as 

P7"-(r, L) = (^^^l-,yi^(^^^^^y {A' exp[-Ai0f (T)] + exp[-A20f (T) + A^H exp[-A3<^3-(r)]} , (62) 
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with the functions <f>f (T) and 0,f (T) given by ([5^ and p3p respectively, and the function (p^ (T) is 

</),^(T) = {L- wi)2 + (L - t;2T)2 + t[{vl - vl)T^ - 2L{vi - V2)T]T/T,^^ - lAa^Ei - E2)T + tW,{pi -p2)L. (64) 

where we have neglected the terms with powers higher than 0{T/T,p^). 

Now we focus our attention on the elimination of the time dependence that is present in the neutrino oscillation 
probabilities. To do it, we take the average on the time of the expressions (|62|) and (|63|) . obtaining 

where the integrals in the time /i , I2 and I3 are 
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(70) 

The integrals /i and /2 can be performed using both Gaussian integration and the Laplace approximation method, 
while the integral /a can be only performed using the Laplace approximation method. After the time integrations are 
performed, we obtain from (|65p and ()66p the following neutrino oscillation probabilities 
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with the functions f^, f^ and /g^ given by 
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In the relativistic limit, using the approximations (|42|) . (|43)) and (|44)) . we obtain from (|7ip and ((72)) that the neutrino 
oscillation probabilities are 
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Pj^"'^(L) = l--sin^[20i2]n 



1 



(1-^4)1/2 



exp 



i27r- 



L 



L' 



L \ _^ 2c2 f'^ri^ - a-a)^ 



coll 



Los 



(79) 



n?"'"(i) = Jsin2[2Mil 



(1-04)1/2 



exp 



i27r- 



coh 



o 2c2 /'^'■(l -as) 



(80) 



where £' ^ and i' ^ are written as 



T' =_ 

osc 



L' 



\- ai 

Lcoh 



coh 



(1+02)1/2 



and 



fli 



0-2 



0-3 



U'(a™?2 



8 E^ 

5^ ^2' 
5 4(Ato22)2^2 



16 



£;6 



a4--- 



3 4(ATOf2)^'^^ 



^6 



(81) 
(82) 

(83) 
(84) 
(85) 
(86) 



We can observe that the oscillation probabilities (|79p and ([50)) have the same functional form that the oscillation 
probabilities (|45|) and (|46l) obtained considering only the space spreading of the wave packets, but now the exponential 
is multiplied by a factor that includes 04. But new, due to the time spreading effects of the wave packets, the 
expressions ([5T|) and ((5^ show respectively some changes of the oscillation length (|T7l) and of the coherence length 
p8)) . We observe in the space-time spreading case how the oscillation length ([8T|) and the coherence length ([82]) are 
respectively a little smaller than the ones obtained in the space spreading case given by (|47|) and (|48|) . Due to the 
time spreading effects, now the maximum number of oscillations is 



A^' 



L' 



coh 



(l + ai) 



T I 

^osc 



(1 + 02) 



AT 

1/2 °«'=' 



(87) 



which implies that it is smaller than the one obtained in the space spreading case (see the expression (j5ip V However, 
working in the regime of oscillations with transversal spreading and for the limits of nearly degenerate masses and 



relativistic neutrinos, the quantities Oi, with i = 1,2,3,4, are very small, so L'^^sc — Lose, L'^^y^ ~ Lcoh, 



Ksc ^ No 



as ~ and 04 ~ 0. In this way, the oscillation probabilities ([75]) and (|M)) can reduce to (|^ and (|i5)) . Thus, we find 
for relativistic neutrinos that the time spreading effects of the wave packets in the neutrino oscillation probabilities 
can be neglected in the regime of oscillation with transversal spreading. 



VI. CONCLUSIONS 

We have studied the space and the time-space spreading effects of the wave packets on the neutrino oscillation 
probabilities in vacuum for considering that neutrinos are described by Majorana fields. We have calculated the 
two-flavor neutrino oscillation probabilities in the context of a wave packet extension of the quantum field theory 
treatment that we previously developed for the case in which the Majorana neutrino mass eigenstate were described 
by plane waves. In the treatment that have presented here, the flavor neutrinos have been considered as superpositions 
of neutrino mass eigenstates described by localized wave packets. We have studied the spreading effects of the mass 
eigenstate wave packets on the two-flavor neutrino oscillation probabilities in the case of considering only the space 
spreading and in the case of considering simultaneously the time and space spreading. 
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We have showed for the space spreading case, in the relativistic hmit, that the Majorana neutrino oscillation 
probabilities have the same functional form as those reported in the literature for Dirac neutrinos. On the other hand, 
for the time-space spreading case, in the relativistic limit, and working in the regime of oscillations with transversal 
spreading, we have obtained a new expression for the two-flavor neutrino oscillation probabilities that includes some 
small changes in the oscillation and coherence lengths. Due to the time spreading effects, the oscillation length, 
the coherence length and the maximum number of oscillations obtained in the space-time spreading case are a little 
smaller than the ones obtained in the space spreading case. These differences can be neglected in the relativistic limit, 
and consequently the neutrino oscillation probabilities of the time-space spreading case can reduce to the ones of the 
space spreading case. 
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